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We discover presence of a hitherto unexplored type of resonance in a parametricaUy excited Van 
der Pol oscillator. The oscillator also possesses a state of anti-resonance. In the weak non-linear 
limit, we explain how to practically get a complete picture of different states of limiting oscillations 
present in the oscillator when the non-linear term therein is excited by an arbitrary 2n periodic 
function of time. We also illustrate how two such oscillators can be coupled to behave like a two- 
state switch allowing an sharp change of value of amplitude for stable oscillations from one constant 
to another. 
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Oscillators — linear and non-linear — are ubiquitous 
in mathematical models in almost any research area and, 
thus, resonance is one of the most common phenomena 
around us. A simple harmonic oscillator (or any oscil- 
lator, for that matter) or a system of such oscillators is 
said to be at resonance when it oscillates with far greater 
amplitude at certain frequencies (resonant frequencies) 
than at others in the presence of an external driving pe- 
riodic force which not necessarily be large. Even in ab- 
sence of any driving force, the simple harmonic oscillator 
can have exponentially growing oscillations under certain 
conditions, for instance, when it is parametrically excited 
meaning that its parameter oscillates with time. Such 
parametric resonance [1] is well studied in the context of 
Mathieu equation appearing across disciplines over the 
last century. Also of interest and far more generality and 
complexity are the nonlinear parametric equations of the 
form [2]: 

£ + k±+[ujl+p{t)]^x + J2a^x^+'^ =0. (1) 

Here, dot denotes derivative w.r.t. time, wq is the natu- 
ral frequency of the oscillator, p{t) is a periodic function 
and k the damping coefficient. Damped Hill equation is 
recovered when ai = 0. While for p{t) (x cos(2cjot) in 
Hill equation, we have the well-known (linearly damped) 
Mathieu equation. These parametric excitations, when 
favourable, set up instabilities in the center-type oscilla- 
tions (i.e., amplitude of the unforced oscillation is differ- 
ent for different initial conditions). 

For exploring the effect of parametric excitation on 
self-sustained limit-cycle-type oscillations, Lienard equa- 
tion given by 

x + f{x)i + gix)^0, (2) 

where f{x) and g{x) are continuously differentiable func- 
tions on real line, is a natural candidate for investigation. 
This is because under suitable conditions, Lienard the- 
orem (see e.g. [3]) guarantees a stable and unique limit 



cycle surrounding the origin in Lienard systems. One of 
the simplest and very frequently studied Lienard system 
is the Van der Pol oscillator [Eq. (3) withp(t) = 1]. Van 
der Pol equation, along with its modified forms, have 
proved to be of practical importance in various physical, 
chemical and biological sciences [4-6] as well as in eco- 
nomics [7] , and of late, even in seismology [8] . Parametric 
excitation on Van der Pol oscillator has been researched 
in the context of Mathieu- Van-der-Pol equation [9-12], 
however the parameter in front of the non-linear term in 
the equation remains constant therein. This constraint, 
as we shall see below, when relaxed opens up a plethora of 
hitherto unexplored exciting possibilities such as an un- 
conventional resonance and anti-resonance phenomena. 
In passing, it may be mentioned that anti- resonance [13], 
occurring in multi-degree-of-freedom systems at frequen- 
cies in between two resonating frequencies, is an impor- 
tant phenomenon studied in the context of seismic pro- 
tection of structures. Also, among many other appli- 
cations of anti-resonance, driving piezoelectric motor at 
anti-resonant frequencies has practical advantages [14]. 

Here, we embark on a study of an equation for the 
"parametrically excited nonlinearity in Van der Pol 
ocsillator" ("penvo" hereafter) defined as 

X + ep{t){x'^ - l)x + uj^x , (3) 

where p{t) is a 27r-periodic function and, hence, possesses 
following Fourier series: 

oo oo 

Pi^) = Y + X! cos{nt) + '^bn sin {nt) . (4) 

n— 1 n—1 

Very interesting results, that form the core of the study 
presented here, and illuminating insights can be gained 
by using a rather simplified truncation of the series (4): 
p{t) — aQ/2 + a„ cos(nt). Below we write down the equa- 
tion for the simplified penvo thus obtained: 

X + e[l + cos{nt)]{x^ ~ l)x + X ^ . (5) 

Here, 7 is a real number. We have rescaled time t and 
parameter e in order to fix wo = 1 and ao = 2 without 
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FIG. 1. Anti-resonance and resonance in simplified penvo: 
= 2 {left) and f2 = 4 (right). The solid line in the left fig- 
ure is given by •\/(4 — 2|7|) and by 1/8/(2 — I7I) in the right 
figure. The dot-dashed and the dashed lines have the constant 
values of ^/2 and 2 respectively. Circular markers which col- 
lectively appear as dense thick lines are the numerical data 
points in excellent fit with theoretical predictions. 

any loss of generality. Also, we have generalised integer 
n to a positive real number f2 in order to increase the 
scope of our study. We shall come back to penvo [Eq. 
(3)] again at a later stage. We remark here that we shall 
solely work in the weak non-linear limit; |e| <C 1 and 
positive e where the limit cycle in Van der Pol equation 
[Eq. (5) with 7 = 0] is stable. Unless otherwise specified, 
e has been fixed at 0.001. 

We begin with briefiy introducing the main numerical 
results that are typical of the simplified penvo. For any 
il / 2 or 4, double precision numerics show that Eq. (5) 
allows for a solution corresponding to a stable limit cy- 
cle of radius 2 — accurate upto 5% of relative deviation. 
When Q is irrational, as one may guess, quasiperiodicity 
can also be seen in the stable oscillations. J7 = 2 gives 
rise to oscillations of amplitude less than 2 and decreasing 
with increasing I7I until I7I is unity, after which ampli- 
tude saturates at a value •\/2(< 2). This can be, aptly, 
considered as an anti-resonance phenomenon. At £7 = 4, 
the simplified penvo shows resonance: as I7I is increased 
from to 2 the oscillations' amplitude increase mono- 
tonically and always remain greater than 2. Readers are 
referred to Fig. 1. 

It is worth noting that while to see the resonance one 
doesn't need an external driving force, this is also unlike 
the parametric resonance seen in the Mathieu oscillator 
in at least following two respects: (i) the resonance in 
simplified penvo is not an instability phenomenon where 
amplitude of the oscillations increases with time, and (ii) 
for resonance the parametric periodic forcing has a period 
one- quarter of the natural period of the simplified penvo 
where the angular frequency of the limiting oscillations 
for all 17 is 1 up to the first subleading order in e. 

In order to explain the numerical results, we resort to 
the harmonic balance method [15]. We assume a highly 
truncated Fourier series expansion of the periodic solu- 
tions under question: 

x{t) = a{t) cos(t) + b{t) sin(<) , (6) 

where we assume slow enough time dependence of co- 
efficients a and b so that we can ignore a, b, ed and eb 
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FIG. 2. Phase plots for simplified penvo {Q = 2) in the Van 
der Pol plane — 7 = 0.5 (left) and 7 — 1.5 (right): The 
fixed point 5* is absent for < 7 < 1 when R is the stable 
fixed point but shows up as the sole stable fixed point in the 
positive quadrant when 7 > 1. 

(|e| << 1) w.r.t. d and b. After substituting Eq. (6) 
in the simplified penvo [Eq. (5)], one equates the co- 
efficients of cos(f) and sin(t) to zero. Consequently, we 
can find the flow equations for a and b. The trajectories 
of the flow equations are, thus, easily visualized in the 
Van der Pol plane i.e. the phase plane for a and b. It is 
convenient to consider three distinct cases. 

First let us focus on the interesting case of = 2 in 
details. The flow equations are 

^^f[('-i)4-(i-i)'i 

There are nine fixed points for Eqs. (7) which are given 
by: (a^b^) = (0,0); (0,4 + 27); (4 - 27,0) and (1 + 
1/7,1 — 1/7). The pair of eigenvalues of the characteristic 
matrix obtained by linearising Eqs. (7) about these fixed 
points respectively are: e(2 ± 7)/4; £7(7 + l)/2, —£(7 + 
2)/2; £7(7-l)/2,£(7-2)/2and -£/4±£(5 - 472)1/2/4. 
Particularly, in the first (positive) quadrant of the a — 6 
phase plane, the real fixed points for I7I < 1 are (cf. 
Fig. 2) three points P, Q and R given by (0,0), 
(0, +(4 + 27)1/2) and (+(4-27)1/2,0) respectively. P 
is an unstable node, Q is a saddle point and R is the 
only stable fixed point. The amplitude of oscillation cor- 
responding to i? is (4 — 27)1/2 which we have used in 
Fig. 1 to fit the numerical data. (In passing, we men- 
tion that when 7 < 0, instead of R, Q becomes the only 
fixed point in the first quadrant; but the correspond- 
ing stable oscillation's amplitude is still (4 — 2|7|)i/2.) 
However, the fitting breaks down after I7I > 1. So, 
what happens at I7I = 1? Assume 7=1 (argument 
for 7 = — 1 goes along similar line). At 7 = 1, another 
real fixed point S given by ((1 + 1/7)^/2, (1 - 1/7)^/2) 
emerges coincided with R. As 7 becomes greater than 
1 , an inspection of the corresponding eigenvalues tells us 
that among P, Q, R and S, only S' is a stable fixed point 
in the positive quadrant. The emergence of new stable 
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fixed point fixes the amplitude of the stable oscillations 
at (a2 + 62)1/2 ^) = (1 + 1/^ + 1 _ 1/^)1/2 = ^ for 

all 7 > 1. 

The case = 4 can be explained similarly. We have 
in the Van der Pol plane trajectories governed by 

o = - [8 - (2 + 37)5' + (7 - 2)a2] (8a) 
^ = § [8 - (2 + 37)a' + (7 - 2)6^] . (8b) 

There are nine fixed points for Eqs. (8) which are 
given by: {a\h^) = (0,0); (0,8/(2 - 7)); (8/(2 - 
7), 0) and (4/(2 + 7), 4/(2 + 7)). The pair of eigenval- 
ues of the characteristic matrix obtained by linearis- 
ing Eqs. (8) about the above fixed points respectively 
are: e/2, e/2; -2e, 476/2(7 - 2); -2e, 476/2(7 - 2); and 
— e, 276/2 -|- 7. With the help of these information, one 
can show that the amplitude of the periodic oscillations 
is 2\/2/{2 — — a function which fits well with the 

numerical data presented in Fig. 1. For I7I > 2, there is 
no stable fixed point of Eqs. (8). 

Finally consider the case: O ^ 2,4. When Vi is an 
integer, the flow equations in the Van der Pol plane are 

« = y(i-T-4j = ^ = y(i-T-4j-(9) 

The origin (a, b) ~ (0, 0) is a fixed point of Eqs. (9). The 
pair of eigenvalues of the characteristic matrix obtained 
by linearizing Eq. (9) about the origin are: e/2,e/2. 
Encircling this unstable origin, there exists a stable limit 
cycle in the Van der Pol plane: a^ + b^ = 4. This explains 
the fixed amplitude of 2 for the limit cycles in the simpli- 
fied penvo when f2 takes integer values other than 2 and 
4. For any non-integral value of O, the harmonic balance 
procedure adopted here conceptually breaks down. How- 
ever, straightforward application of any standard pertur- 
bation technique such as renormalization group method 
[16, 17], Lindstedt-Poincare method (used below) etc. 
tells us that the corresponding limit cycle's radius re- 
mains 2. 

Having completely explained the amplitude of the sta- 
ble oscillations in the simplified penvo, it is straightfor- 
ward to find out the correction to the natural frequency 
{ujQ — 1) using Lindstedt-Poincare method [15]. To this 
end, we assume following perturbation series: 

x{t) = xo{t) + exi{t) + e^x2{t) + 0{e^) , (10) 
Lu{t) = LOo{t) + euji(t) + e^uj2{t) + 0{e^) . (11) 

Usage of these series in Eq. (5) leads to following solution 
for xo{t): 

xo = A cos{ujt) + B s\x\{ujt) (12) 

using which in the next order equation for xi{t), and 
subsequently equating the coefficient of cos{t) and sin(t) 
to zero in order to get rid of secular terms, we arrive at 
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FIG. 3. Quasiperiodicity in simplified penvo (fi — -\/2): {left) 
A trajectory is seen to densely fill up a portion of x — i plane 
of phase space, and {right) an elliptical curve traced by the 
same trajectory in the relevant Poincare section (cf. text). 

two simultaneous conditions on as a function of A and 
B. One can easily argue that — for the periodic oscil- 
lations found using harmonic balance method — uii ~ Q 
for all VL. Thus, the frequency of the stable oscillations 
in simplified penvo, to the first subleading order in e is 

w = l+0(e2). (13) 

With a view to exploring the existence of quasiperiod- 
icity, let us rewrite Eq. (5) as an autonomous set of three 
first order ordinary differential equations: 

±1= X2, (14a) 
±2 = - £(1 + cosa;3)(l - x2)2;2 , (14b) 

±3=0. (14c) 

Whether for irrational J7, the oscillations are quasiperi- 
odic or not, is best addressed by looking at the Poincare 
section at some plane: X3 =constant. Since the third 
phase space coordinate X3 is the phase of the parametric 
excitation, the appropriate Poincare section corresponds 
to looking at the oscillations after every period 2ti/Q,. 
For periodic oscillation, depending on the ratio between 
angular frequency of the oscillation and O, the Poincare 
section will consist of finite number of discrete points, 
whereas any quasiperiodic dynamics will generate infi- 
nite number of points on the Poincare section so that a 
dense curve will be seen — this exactly is what we see in 
simplified penvo with irrational O (cf. Fig. 3). 

Let us come back to Eq. (3) with ujq = 1 without 
any loss of generality. First consider the case of p{t) 
being an odd function so that the Fourier series (4) of 
p(t) is an exclusively sine series. One can easily show 
using harmonic balance — as used before — that for such 
a form of p{t), there doesn't exist a single stable node 
or spiral in the Van der Pol plane. Consequently, the 
limit cycle in Van der Pol equation is destroyed and there 
is no limiting periodic oscillation. When p{t) is a even 
function, all the sine terms in the Fourier series of p{t) 
vanish. Rescaling e and defining ^ci = 2ai/ao, we find 
that, to the leading order, only 7c2 and 7c4 can determine 
the conditions for stable periodic oscillations. Other 7ci's 
are irrelevant for our purpose: when 7c2 = 7c4 = 0, penvo 
will always yield a limit cycle solution — + b^ = A — 
irrespective of what values other 7ci's take. 
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In the most general case, a complete description of 
penvo in the oscillating state effectively requires any ar- 
bitrary 27r periodic p{t) to be represented by the following 
surprisingly small fraction of the full Fourier series (4): 

ep{t) = £ (1 + 7c2 cos 2t + 7c4 cos 4t + js2 sin 2t + 7s4 sin At 
-|-irrelevant terms) , (15) 

where ao/2 has been absorbed in e and 7si = 2bi/aQ. 
When only 'l+irrelevant terms' in Eq. (15) survive, the 
usual Van der Pol oscillator's limit cycle is recovered. 
While dealing with all four 7's simultaneously is cum- 
bersome, we can understand relative importance of the 
sinusoidal terms in expression (15) by setting any two of 
the four 7's to zero at a time. Out of the six cases that 
can thus arise, only the following forms oip{t) require in- 
dependent study: (i) 1 + 7c2 cos(2t) -|~ 7c4 cos(4i) and (ii) 
1 + 7s2 sin(2t) 4-7s4 sin(4t). Remaining four relevant com- 
binations of sine and cosine can be reduced to these two 
forms by suitable rescaling of time t and the concerned 
7's. Since we do not find any qualitatively new result 
other than already presented earlier, we merely mention 
that though the presence of four parameters in p{t) make 
algebra tedious, one can in principle go through it and 
theoretically explain the numerical results in a straight- 
forward fashion along a similar line of argument adopted 
herein for the simpler cases. 

Now we present an idea how a two-state-switch may 
be constructed using a coupled system of two simplified 
pcnvos. Consider, 
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FIG. 4. Two-state switch in coupled simplified penvos. The 
behaviour of switch for j7| > 1 is very sensitive to fluctuations 
in the value of '2y/2/ Ay^' . To compensate this effect, Ayoo 
has been calculated by taking average of 100 points over the 
limiting cycles of corresponding oscillations. The dashed line 
is for amplitude^ 2, the dot-dashed line is for amplitude= y/2 
and circles are the numerical data points. 
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When < I7I < 1, Eq. (I6b) results in an oscillation of 
amplitude Ay^ = {yl^+yloY^^ = (4- 2|7|)i/2 such that 



Ay^ g [2, V2), where subscript 00 refers to the quantities 
calculated at i — > 00. This in turn will mean that the 
angular frequency in the cosine term in Eq. (16a) lies 



between y/2 to 2 (but not 2 yet). Hence the amplitude 
of the stable oscillation one obtains from Eq. (16a), for 
7 € [0,1), is always 2. When I7I > 1, Eq. (16b) yields 
a periodic (or quasiperiodic) oscillation of amplitude \/2 
and hence the angular frequency in the cosine term in 
Eq. (16a) is fixed at 2. Consequently, the amplitude of 
the stable oscillations one obtains from Eq. (16a) is fixed 
at for I7I > 1. In short, the amplitude of the stable 
oscillations generated by Eq. (16a) switches from 2 (for 
I7I < 1) to (for I7I > 1) at I7I = 1 (cf. Fig. 4). 

To conclude, we emphasize that in view of this study 
on Van der Pol oscillator, it is of immediate interest to 
revisit different Lienard systems displaying limit cycle 
oscillations and investigate the effect of parametric exci- 
tation on the non-linear terms. We suspect that results 
obtained here could be generic and, the resonance and 
the anti-resonance could be encountered in other such 
systems as well. Such resonance and anti-resonance can 
have direct implications on every system across disci- 
plines that uses one or the other form of such oscillators 
for its mathematical modeling. 
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